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Abstract 

Some new identities for quantum variance and covariance involv- 
ing commutators are presented, in which the density matrix and the 
operators are treated symmetrically. A measure of entanglement is 
proposed for bipartite systems, based upon covariance. This works 
for two- and three- component systems but produces ambiguities for 
multicomponent systems of composite diemnsion. Its relationship to 
angular momentum dispersion for symmetric spin states is described. 

1 Introduction 

Several measures of entanglement |l| or quantum correlations have been pro- 
posed: some are associated with the preparation of the state, others with 
the process of purification or distillation [Q and yet others with the notion of 
mutual information or relative entropy In this paper we wish to suggest 
another measure, based on covariance, in which the acts of state creation and 
observation are considered in a dual manner. 

In practice it is very natural to describe the condition of the system (its 
method of preparation or lack of it) in terms of a density matrix which is 
tied to the subsequent observations on it. This is how the linkage between 
observer and observed occurs quantum mechanically, and of course the results 
are expressed in terms of traces over appropriate functions of the density 
matrix and of the operators being measured |Q. Indeed Mermin ^ has taken 
that view that the density matrix, and the correlations between observables 
which thereby ensue, constitute all of physical reality. 
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In this paper we will also focus on the density matrix. Because binning 
of observations is a necessity in practice, the dimension of the density matrix 
is thereby determined: N separate bins produce a density matrix p that is 
an X hermitian matrix, satisfying the usual hermiticity and trace condi- 
tions. In this way, we can regard the basis as an A^-level system, rather like 
a particle of angular momentum J = {N — l)/2. Hence, although we might 
be studying the probability distribution of an observable which actually pos- 
sesses a continuous spectrum, we can still regard it as a spin-like system; in 
practical terms, the bigger the binning number N, the greater the precision 
of the information about the continuous variable, but obviously is never 
infinite. For spin measurements, we need not go to such pains because N is 
fixed for us at the start. 

With the focus on density matrices, we will carry out measurements (with- 
out mutual interference) on two subsystems, 1 and 2 say, so their correspond- 
ing observables, superscripted by (1) and (2), are commuting operators. The 
system will be separable ^ or "disentangled" if the larger density matrix p is 
merely a direct product of density matrices associated with the two subsys- 
tems, or p = p*^^) ® p*-^-* ; a particular case arises when the initial or prepared 
state is the direct product of two subsystem states, {ip) = \4'^^^)\x^'^^)- When 
two subsystems are disentangled, the results of measuring any quantity 
in the first subsystem are not tied to the results of measuring any quantity 
B^"^^ in the second subsystem; necessarily 

for all choices of A and B. However, if p 7^ p*-^' ® p^^-*, the configuration is 
non-factorizable and the covariance, 

cov(A(i)5(2)) = - (1) 

no longer disappears. 

The real issue is how to quantify the entanglement or lack of factorizability 
[^] of the larger density matrix. Several proposals have been advanced in the 
literature, but none of them is entirely simple or definitive However all 
researchers in this field seem to agree on the following three conditions for 
an entanglement measure E{p): 

1. E{p) = iff p is separable, ie if the density matrix can be written as 

(1) ^ (2) 
P = 2^iP^P^ ® Pi ■ 
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2. Local unitary transformations should leave E{p) invariant. 

3. E{p) should not increase under local measurement and classical com- 
munication procedures, we intuitively know that such procedures can- 
not add non-locality characteristics to the system being measured. 

As an extra requirement, it would be nice if E{p) gave some indication of 
the extent of violation of Bell-type inequalities |^ . 

In this paper we want to put forward a concrete scheme for quantifying 
correlations between two subsystems and their possible entanglement. The 
scheme is based on a generalization of eq. (1) and particular choices of 
operators A^^^ and B^'^\ which are readily applicable and rooted in the density 
matrix notion. In the next section we discuss several matters connected 
with non-separability of states and their influence on subsequent subsystem 
measurements. Because we deal with practical observations, the density 
matrix is truly discrete and we can assume that the elements of the vector 
space on which it lives have equal weight. As already mentioned, one may 
regard the dimension = A^^^^ -|- A^*^^-* as corresponding to a "spin system" , 
with each component carrying equal weight, and can adopt the same stance 
for the subsystem dimensions A^*-^'^-*. (This restriction can be relaxed if the 
components have unequal weights, such as atomic energy levels at a finite 
temperature.) 

The next section deals with the generalities of simultaneous measurements 
and their covariance properties |T^. This is followed by our suggestion for 
quantifying entanglement of two subsystems within a larger entity, which is 
shown to be consistent with normal expectations for two spin 1/2 subsystems, 
when A^^^) = A^*^^^ = 2. We also discuss the use of total spin dispersion [0 
as another measure of entanglement, with an allied appendix concerning the 
Majorana- Penrose ^ representation of spin states on the Poincare sphere. 
The subsequent sections deal with entanglement measures for larger value of 
A^^^) and N^'^\ Finally we discuss general questions pertaining to our sug- 
gested measure; these include Rovelli's notion that information in quantum 
mechanics is relational [Q, Mermin's notions of correlations between local 
observables , and the difference between our modified correlation measure 
with classical correlations for impure states. 
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2 Correlations and density matrices 



Elementary texts on quantum mechanics teach us that the results of all 
physical measurements and processes can be tied to the evaluation of traces of 
products of observables with the hermitian density matrix p. Thus statistical 
formulae like 

(F) = Tr(pF); Tr(p2) < 1, 

etc. are part of the standard repertoire. Of course, the p-eigenvalues lie 
between and 1; in the latter case we are dealing with a pure state when 
the density matrix reduces to a projector p ^ = while the most 

random situation p = 1/N corresponds to the case of maximum entropy. 

The covariance for any two commuting observables A,B in a mixed state 
p is defined as 

coVp(A, B) = {AB) - {A) {B) = tiipAB) - tr(pA)tr(p5). (2) 

Clearly, varp(y4) = cov p{A, A). Less well-known is the fact that pure state 
dispersions and correlations can be neatly expressed in terms of a single trace. 
Consider the quantity 

B) = tr([p, A][B, p])/2 = tr(p2{A, B}/2 - pApB), (3) 

where A and B are any two operators. This quantity will be referred to as 
the alternative covariance Q 

We now present some elementary results about Cp which follow simply 
from this definition: 

1. Cp{A, B) = Cp{A — a, B — b) , where a, b are any two constants. 

2. Cp{aA,bB) = abCp{A,B). 

3. Cp{E^^A,Z,B,) = Z^,Cp{A-Bj). 

4. Cp{A,A) = tr(pM2 - pApA). 

^Evaluating traces of larger numbers of pure state commutators, one may establish 
algebraically that for odd numbers of products, the traces do vanish. For instance, 
tr([Ap][S,p][C,p]) = 0, etc. 
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5. C^p^t{A, B) = Cp{U'^AU, WBU). where U is any unitary transforma- 
tion. Thus a change of basis for the state is equivalent to an inverse 
change of basis for the operators. 

6. Cp{A,A^)Cp{B,B^) > \Cp{A,B^)\\ This fohows by considering the 
operator T ^[p,A- cB], with c = tr([p, A][A\ p])/tr([p, A][B\ p], and 
noting that tr(TTt) > 0. 

7. Cp{A, A^) is real. 

8. \Cp{A, B)\'^ is symmetrical under interchange, conjugation and change 
of phase of the two operators. 

All of these properties are shared by the usual covariance cov p{A, B). Never- 
theless, alternative covariance Cp does not provide an indication of variance 
and covariance in the usual sense. For instance, if the state is one of maxi- 
mum entropy, on the one hand we have Cp{A, B) = for all A, B because p 
is proportional to unity; on the other hand, the coVp{A, B) need not vanish. 
Some special cases for the operators A, B can now be studied. 

1. If A and B commute, Cp{A, B) = tr(pMS - pApB). 

2. If A and B are both hermitian, Cp{A, B) becomes real. 

3. If A and B arc both unitary, Cp{A, A^) = tr(p2 _ pApA^) < tr(p2) < 1. 
Likewise for B. Since Cp{A, A^)Cp{B, B^) > \Cp{A,B^)\\ it follows 
that \Cp{A,B)\'^ < 1. 

More particular cases arise when the system is prepared in a pure state 
so that p becomes a projection operator and Cp{A, B) reduces to 

Cp{AB) i(v;|{A5}|^)-(^l^l^)(^|5|^) = ^({^-(^),i?-(i?)}) 

= cow p{A,B), when [A,B] = Q. (4) 

Thus 

Cp{A,A) ^ {A') - {Af ^vs..p{A). (5) 
This is in keeping with the familiar variance-covariance inequality: 

var,(yl)varp(5) > \coy p{A,B)\\ (6) 
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If A and B are commuting unitary operators and because \Cp{A, B)\'^ < 
Cp{A, A)Cp{B, B) < (tr(p-^))^ < 1, we see that alternative covariance only 
attains a value of 1 for pure states. 

It is worthwhile comparing the two covariance functions, in relation to 
two commuting observables, A, B. Since [A, B] = 0, select an orthonormal 
basis \i) wherein the operators are simultaneously diagonalised, so 

^ XI ^ ^ \i)bi{i\- 

i i 

Then 

C0Yp{A, B) = ii{pAB) - tr(pA)tr(p5) 

= ^aihi{i\p\i) - ^aihj{i\p\i){3\p\3) 

= ^ai{bi-bj){i\p\i){j\p\j), 

and similarly 

Cp{A,B) = ii{p^AB) -ti{pApB) 

= ^ai{hi-hj){i\p\j){j\p\i), 

since J2j{j\p\j) = tr(p) = 1. Furthermore note that Xlj j(^|p|j)(j|pK) = 
tr(p^) < 1. Upon symmetrising the sums, we obtain the neater expressions, 

coyp{A,B) = ^(a, - aj){h - bj){i\p\i){j\p\j) /2 (7) 

Cp{A,B)^J2(^,-b,){b,-b,){i\p\j){j\p\i)/2 (8) 

Whilst the ordinary covariance has a clear meaning — namely, a measure of 
the correlations between the results of local measurements A and B that 
commute — the interpretation of the alternative covariance is less obvious. 

We can obtain more insight by choosing A = B. Since p is a positive 
definite hermitian operator, {i\p\'i){j\p\j) > (^|pb)(j|pK)) for two states 
Therefore for a general (mixed state) density matrix, 

^{tti - aj)'^{i\p\i){j\p\j) > ^{tti - ajf{i\p\j){j\p\i), 
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or 



M^ip{A)>Cp{A,A). 



(9) 



In the light of the variance inequahty above it is natural to ask whether 



is true. In fact a single (but carefully chosen) counterexample suffices to 
show that it is false: in local bases u, d for two local operators A®1 and 
1 (8) take 

p= \{uu + dd)){{uu + dd)\/A+ \ud){ud\/A+ \du){du\/A 

and select the local operators to be diagonal, 

A®1 = {\uu) {uu\ + \ud) {ud\) — {\du) {du\ + \dd) {dd\) 
1 ® i? = {\uu) {uu\ — \ud) {ud\) — {\du) {du\ — \dd) {dd\) . 

Evaluation of the two types of covariance leads to 



Thus the variance inequality cannot be extended to covariance. 

However, an immediate consequence of the inequality, varp(A) > Cp{A, A), 
is that when vaip{A) = 0, Cp{A, A) = too for any observable A. But 
2Cp{A,A) = ti{[A, p][A, pY); so [A, p] = 0, which means that p is purely in 
an eigenstate of A. This accords with the basic tenets of quantum mechanics 
of course. The contrapositive of this result is that if [^4, p] 7^ 0, then varp(A) > 



Another worthwhile comment stems from the observation that if X is 
conjugate to A in the sense [A, X] = ih, then 



COY p{A,B)\ > \Cp{A,B)\ 



coy p{A, B) = 0, but Cp{A, B) = 1/4. 



0. 



2Cp{A,A) = tr{[p,A][p,A]^)^hhr 



( 



dp dp ^ 
dXdX 



) 



Thus, 
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with equality only applying to pure states. For example, the energy uncer- 



|2 



tainty is given by (AH) > h tr (^|^| j /2, while the momentum uncertainty 
is given by the derivative of the density matrix with respect to position: 
{APy > ^2 tr (^lllf) /2, and so on. 

For a general mixed configuration, the two inequalities, 

varp(A)var^(5) > Cp{A, A)C,{B, B) > \C,{A,B)f 

together with the well-known 

YaTp{A)vaTp{B) > \cov p{A, B)\'^, 

provide a lower bound for the experimentally observed variance products of 
any two operators, whether or not they commute. For instance, 

varp(X).varp(P) > Ihhi ^ 



dXdP 

In the next section we present examples of operators A, B, for which there 
exist states such that \covp{A, B)\ > \Cp{A,B)\ and also other states for 
which \Cp{A,B)\ > \coVp{A, B)\. Hence &ot/i inequalities must be considered 
jointly in an examination of the minimum of the variance product, together 
with Heisenberg's well-known lower bound, |tr(p[A, S])|/2. 



3 Correlation measures for pure states of two 
subsystems 

This section examines the correlation properties of the entanglement of two 
subsystems in a tensor product Hilbert space H(i) ® H(2). By definition, 
measurements can be carried out without mutual interference on the two 
subsystems so their corresponding observables, superscripted by (1) and (2), 
are commuting operators. As mentioned in the introduction, a state of the 
system is factorisable or disentangled if the larger density matrix p is merely 
a direct product of density matrices associated with the two subsystems, or 
p = p^^^ (8) p^^^. For any two local operators A^^^ — A^l, B^^^ = 1 5, it is 
easy to show that the covariance disappears: 

cov,(i)^p(^)(^^'\ B'^'^) = {{A - (A)) ® {{B - {B))) = 0. (10) 
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This includes the case of a pure disentangled state, |0) = \(f)'^^)\(p^'^). 

Having noted that the covariance is non-zero in disentangled states, we 
now refer to the conditions imposed upon any measure of entanglement. The 
second condition is that it be invariant under local unitary transformations. 
With this in mind, define the covariance entanglement for pure states as 

= max \covupwiA^'^\B^'^^)\. 

The maximum will clearly be invariant under additional local unitary trans- 
formations. 

Since the operation of permuting the elements of a Hilbert space is uni- 
tary, all elements of the Hilbert spaces are equally important. For this reason, 
it is natural to select the operators A^^^ and B^"^^ so as to equally weight the 
elements of the Hilbert space. The next section describes several ideas for 
achieving this, starting with the simplest case. 



3.1 Pure state correlations for iV(i) = Ar(2) = 2 



This section investigates a method for quantifying pure state entanglement 
in the simplest possible non-trivial case, corresponding to two spin 1/2 sys- 
tems, with Hilbert space HI (8) H, where H = is a local Hilbert space, 
with orthonormal basis \u), \d). Consider two local operators which dis- 
tinguish between elements of the local Hilbert spaces. With the aim of 
weighting local basis elements equally, define operators in the product basis 
\uu), \ud), \du), \dd) by 



/ 1 



/ 1 



and so 



^(1)^(2) 



.(1)^(2) 



/ 1 



V 



-1 



-1 



Next consider the pure (normalized but arbitrary) state, 

10) = a\uu) + (3\ud) + j\du) + 6\dd). 
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Figure 1: Equal- weight covariance of the Bell-state \uu + dd)/V2 under 
equal local unitary transformations about two axes. 

Working out coVp{A, B) in this state, it is straightforward to show that the 
covariance is maximised provided that \a\ = \5\ = 1/^/2, /3 = 7 = 0, or 
\P\ = I7I = 1/a/2, a = 6 = 0. Thus one may take the four independent Bell 
states, 

|10) = [\ud) + \du)]/V2, |00) = [\ud) - \du)]/V2, 

|1+) = [\uu) + \dd)]/V2, = [\uu) -\dd)]/V2, 

as the ones that have the largest covariance. (These pure states are also 
known to be the most entangled ones.) Of course they are all local unitary 
transforms of just one of them, say the Bell state, |1+) = \{uu + dd))/\/2, 
with a corresponding p = \uu + dd){uu + dd\/2. If one rotates the operators 
A, B together about the "y-axis" by the same amount, we can get a good idea 
of how the covariance varies with rotation angle; maximization is attained 
when the angle is nil. See Figure ^ 

If a pure state is disentangled, then there is an orientation of the local 
A^^^ which has zero variance. To see this, rotate the state by local unitary 
transformations until the reduced density matrix for the local operator in 
question is diagonalised. Since the initial state was pure and disentangled, 
then it may be represented by a separable projector p = p^^^p^'^\ in which 
both reduced density matrices are projectors. Thus the main diagonals of 
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Figure 2: Equal-weight variance of a parametrised set of pure states, under 
different relative orientations of the state and the local operators. The max- 
imally entangled states have x = 7r/4, and are the only states to attain a 
covariance of 1, whilst the disentangled states (x = 0, 7r/2) are the only ones 
which attain a covariance of zero, for particular relative transformations of 
the local operator. 

p(^^ and p*-^-* can be reduced to a single 1, with Os elsewhere. Hence for 
diagonalised local operators A and B in this basis, 

{A') = {Ar, {B') = {B)\ 

so the variances vanishes. Figure ^ illustrates the behaviour of the variance 
in a 2-variable parametrisation. The horizontal axis variable x parametrises 
a set of pure states which range from disentangled to a maximally entangled 
Bell state, and back to disentangled again, i.e. p = \iIj{x)){iIj{x)\ where 
\ip{x)) = cos(a;)|'U'u) + sm{x) \dd) . The second variable y parametrises ?/-axis 
rotations of the local spin basis (1) associated with A alone. 

These results may be applied to actual spin measurements. If one knows 
that a state is pure, but is not certain of the degree of entanglement, local 
spin measurements can be made in a variety of directions. If the variance 
and covariance of these measurements vanishes in a pure state, then the state 
must be disentangled. 
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3.2 Pure state correlations for N^^^ = iV^^) = > 2 

Many different choices of local operators A^^\ B'^'^^ are possible, and different 
choices will lead to different behaviour of the covariance. However before 
considering two sensible choices, let us note that for A^*^^^ = A^*^^) = 2, one of 
the maximally entangled states can be taken to be the state of total spin |10), 
while minimally entangled states are |11), |1 — 1). Now for these state combi- 
nations maximal entanglement happens to equate with maximal dispersion 
(AJ)^ and zero entanglement equates with minimal dispersion (AJ)^, where 
J stands for total spin. This suggests that for higher spin, some maximally 
entangled states might be found by minimising the total angular momentum 
dispersion and vice versa. This approach towards quantifying entanglement 
is quite interesting in its own right and is pursued in Appendix A, where we 
also tie it to the Majorana-Penrose pictorial view of spin. In Appendix B, 
by contrast, we classify any measure of entanglement via an integrity basis 
for density matrix invariants. 

3.2.1 Pair-discrimination 

Consider two-system Hilbert spaces where the local spaces may each have 
dimension greater than 2. Select local operators A^^^ = A0 1,B^'^^ = 1 ® 
B, where A, B may be expressed in their respective local bases as unitary 
transformations of the following diagonal matrix, 

/ 1 
0-10 


V 

Next, maximise the covariance over all such unitarily-transformed matrices; 
the result is perforce invariant under local unitary transformations of the 
state. Labelling the local bases |aj) and \bi) respectively, where i runs from 
1 to A^, operators like the above discriminate between pairs of elements in a 
local subspace of the full Hilbert space, and treat terms of the form lafik + 
ajbi) as the basic element of entanglement^. 

^Another possibility is to replace all the zeroes along the main diagonal of A or B 
with 1 or —1; this makes the operator U unitary, which means that its variance is simply 



(11) 
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3.2.2 Equal weight unitary operators 

Since we wish to handle all the subsystem states democratically, let us define 
an equal weight local unitary operator as consisting of some unitary trans- 
formation of a diagonal matrix comprising the iV-th roots of unity. (Note 
that these matrices are not hermitian when N > 3 and cannot correspond to 
observables.) Here the local weight unitary matrices in their corresponding 
diagonalising basis, up to an overall phase, are given by 

/ exp(2i7r/A^) \ 
exp(4i7r/A'") 

V 1/ 

Only for the spin ^ x spin | case, do these operators correspond exactly to 
the pairwise local unitary operators used previously. The next nontrivial case 
is iV = 3, or spin 1 x spin 1. In this case one can see that all states which are 
local unitary transformations of the pure state |ai6i -t- 0262 + 03^3) / have a 
maximal covariance of 1. To understand why, rotate the equal- weight unitary 
operators so that the local states ai, 02, and &i, 62, &3 produce eigenvalues 
which are respectively conjugate pairs. This yields (A^^^B^^^) — 1. However, 
(A^^^) = {B^'^^) = 0, so the maximised covariance is 1. 

What of other states, such as |ai6i -|- 02^2) when iV = 3? The following 
theorem gives a necessary and sufficient condition for a state to exhibit Cp — 
1, with respect to these equal weight operators and is in agreement with all 
other pure state entanglement measures. 

Theorem: The only pure states which attain the maximised covariance of 1 
under equal-weight local operators are states which are local unitary transfor- 
mations of |ai&i -|- a2&2 -l- . . . -I- aNh^) / . Any other states exhibit smaller 
correlations. 

Proof: In a tensor product space H*^^^ ® M^^^ , consider the state 

i,j 

where pij = |cijp > and 'YliijPij = ^'^'^ the orthonormal bases |ai), \hj) 
are a complete set of eigenstates for the diagonahsed equal-weight unitary 

1 ^ However the distribution of the ±1 eigenvalues is not self-evident, except for 

the spin spin^ case. 
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operators A^^\ B'^'^\ with eigenvalues = cxp{2ink/N),bj = exjp{2i7ij/N) 
respectively. Recalling the result for such A^^\ B^'^'' that 

we see that a maximised covariance of 1 is only attainable in a state where 
(A) — (B) — 0, whereupon the covariance reduces to 



C4A^'\B('^) 



In what cases is this expression maximised, subject to the condition that the 
mean values of the operators remain zero? 

We are seeking to maximise | ^PijCiibj\ subject to 

^Pij = 1' I ^Pij(ii\ = 0, I '^Pijbjl = 0. 
By the triangle inequality, 

I ^pijaibjl <pii|ai6i| + | ^pijaibj] < ... < '^Pij\aibj\, 

where equality holds at every stage only if all the complex numbers Oj, bi 
have equal and opposite phase. This means that we are pairing Oj, bj such 
that non-zero pij (only for i = j) are associated in a one-to-one manner with 
Ojfoj = 1 for all such pairs, otherwise the parallelism of the complex numbers 
will be lost. This will ensure that 

^PiilaAl = ^Pii = 1- 

i i 

At the same time we have to guarantee that the average values of A and B 
vanish or | — 0; I X^Pij^il — 0- Since X]?^? — Si^* = 0, a sufficient 

condition for this is that for all such pairs, the weightings are equal or pu = 
1 / ^/N ; in other words every and its corresponding bi = a* only occurs at 
most once in the terms with equal non-zero weighting. (Actually one may 
introduce an arbitrary phase into cu without affecting this conclusion, but 
we have chosen not to do so.) 
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Having established that the states which maximise the covariance can 
take the form \aial + 0203 + ...), we should point out that it is not necessary 
for all the terms to be paired up. Consider the case = 4, or spin 3/2(8) spin 
3/2 systems, with bases ai, . . . , 04 and 61, . . . , 64 respectively. It is possible 
to attain maximised covariances of 1 under the equal-weight measures both 
for \^) = \aibi + 1 ^J2 and for |0) = |ai6i + 02^2 + 0,363 + a/p/[) 12. This 
is achieved by choosing eigenvalues so that (A*-^-*) = {B^"^^) = and yet pick 
eigenvalues such that the values of A'^^^B^'^^ in the states |aifei) . . . 10464) are 
all 1. As we are in a 4-dimensional space, wc can achieve this by taking the 
eigenvalue sets {1, —1} and {l,i, —1, —i} for both operators on both states 
lip) and 10) respectively. 

This observation means that the 'equal- weight' C-based measures for N > 
2 are not measures of entanglement, under the standard criteria. Information- 
based entanglement measures, such as the relative entropy, specify that lip) 
is less entangled than |0). It appears that the covariance-based measures of 
entanglement are most useful when dealing with spin spin| systems, since 
in this case there is no ambiguity as to the choice of eigenvalues. Similarly 
for operators of prime dimension, such ambiguities are absent, because there 
is only one way to arrive at a maximally correlated state: the non-uniqueness 
only pertains to composite-dimensional local spaces. 

The above result demonstrates that for pure states, the reduced density 
matrices must be diagonalised in order to maximise the covariance of the 
diagonalised local unitary operators. However, for mixed states it is not at 
all obvious that the reduced density matrices must be diagonalised in order 
to maximise the unitary matrix covariance. This issue will will be examined 
in section 4. 

3.3 Pure state correlations for A^*^^^ 7^ N'^'^^ 

For spaces of differing dimension, such as spin Ix spin 1/2, the covariance 
of the pairwise operators behave just as in the other cases. However, if the 
equal-weight local operators are used, covariances of 1 are not attainable. 
This reflects the fact that the bases have different sizes, and so there is no 
way to pair up elements between the bases in a one-to-one manner so as to 
produce a set of product eigenvalues with the same phase. 

The simplest example which exhibits this effect is a spin ^ (g) spin 1 space. 
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with equal-weight operators 

A^^^ = diag(l, -1) ® 1, B^^^ = 1 ® diag(l, exp(2i7r/3), exp(4i7r/3)) (12) 

As with the proof that the states of maximum covariance are |ai6i + . . . + 
ttnbn), a covariance of 1 is only attainable if {A^^^B^"^^) = 1; we also need 
the complex numbers Pijttihj to have the same phase (where the state has 
Schmidt decomposition |0) = \/¥ij\(^ibj) in the basis of eigenvalues of the 
operators A^^\ B^'^^). Following similar arguments to those used in the proof, 
it is not possible to choose pij so that the direction condition is satisfied; this 
is because no repetitions of Oj or bj values can occur in the set of non-zero 
Pij, since the resulting Pijaibj would not have the same phase. But it is not 
possible to partially pair up the given set of eigenvalues so that the directions 
of the products are the same, by straightforward enumeration of the cases. 
Thus states in this basis cannot attain covariances of 1. 

Clearly many other local operators can be defined which provide a variety 
of different correlation measures for two-system states but none stands out. 

4 Correlation measures for mixed states 

Making a distinction between quantum and classical correlations has proved 
a thorny problem in the study of quantum entanglement. The nature of 
the problem may be seen when comparing the two states, one pure and one 
mixed, which possess the same covariance for ^4*^^) = B^'^^ = a^: 

1 11 

-\uu + dd) {uu + dd\ and -\uu){uu\ + -\dd){dd\. 

The first state is a Bell state and is maximally entangled, whilst the second 
state is a mixture of disentangled projectors, and is normally regarded as 

being disentangled. As both states exhibit correlations, it is natural to ask 
whether the alternative covariance introduced in section 2 provides a way of 
distinguishing between classical and quantum correlations. 

4.1 Distinction between Cp and cov^ 

Take any two commuting local measurements, A^^'' = A®\, B^"^^ = 1 <^ B, 
and define the function 

Eab{p) ^ maxt;|C^,^t(AW,£(2))|2, (13) 
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where the maximum is now taken over all local unitary transformations of 
the general mixed density matrix p. In previous sections we investigated the 
behaviour of this function for pure states (when Cp reduces to the covariance), 
and found that it appeared to have many of the properties desirable in a pure 
state entanglement measure. 

The situation where the density matrix corresponds to an impure config- 
uration is more intriguing. Shown below is a comparison of the behaviour 
of the two maximised covariance entanglement measures, in several exam- 
ple configurations, which illustrate the distinction between |coVp 
and |Cp(A«,5(2))|2. 



Density matrix 


|covp(A(i),5(2))|2 




pi = [\uu){uu\ + \dd){dd\]/2 


1 





P2 = [^\uu) {uu\ + ^\uu + dd){uu + dd 


1] 3/4 


1/4 


P3 = \uu){uu\ 








P4 = \uu + dd){uu + dd\/2 


1 


1 



The illustrative state p^ is pure but entangled, ps is factorizable and there- 
fore disentangled, while the matrix pi is not factorizable but can be expressed 
as a sum of separable projectors; therefore pi should represent a disentangled 
configuration, according to standard expectations. By inspecting the table 
we see that coVp{A^^\ -B*-^-*), being nonzero, is not a good entanglement mea- 
sure -E(p), but the alternative Cp{A^^\ B^'^^) is better in that it does vanish. 

The alternative covariance C in mixed states is bounded above by tr(p^), 
so the state must be pure to obtain a covariance of 1 under local unitary 
transformations. Another point worth remembering is that for configurations 
of maximum entropy where p oc 1 and commutes with all operators, C is 
automatically zero. 

Figures ^ and || depict the squares of the equal weight classical covariance 
and alternative covariance of a parametrised mixture of different Bell states. 
They show that both covariance measures attain a maximum of 1 only for 
the pure states, and that the behaviour of these functions depends upon the 
kind of Bell-states involved in the mixtures. 

For any choice of A^^\ B^'^\ a maximised Cp{A'^^\ B^"^^) of zero occurs for 
many mixed configuration p which are disentangled, according to the stan- 
dard definition of a mixed state, such as the states of maximum entropy with 
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(i) |cov,(^W,S(2))| 

i| 

0.8 
0.6 
0.4 
0.2 



0.2 0.4 0.6 0.8 1 

(ii) 

1 

. 8 
0.6 
0.4 
0.2 



0.2 0.4 0.6 0.8 1 

Figure 3: (i) Equal-weight covariance magnitude of a parametrised mixture 
of different Bell-states, (ii) Equal weight alternative covariance magnitude 
for the Bell-state mixture. The paramctrisation in both cases is any one of 
the following: p = a;|6i)(6i| + (l — a;)|62) (^2!, where either = \ud:Ldii) / \/2^ 
I&2) = \ud^du)/^/2, or else = \uu±dd) /\/2, I62) = \ud^ du) /\/2. The 
Bell-states must be different in these parametrisations, to obtain a non-trivial 
result. The result is not a maximisation over all local unitary transforma- 
tions. 
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(i) |covp(>lW,S(2))| 



1 




0.2 0.4 0.6 0.8 1 



(ii) 



1 




0.2 0.4 0.6 0.8 1 



Figure 4: (i) Equal- weight covariance magnitude of a parametrised mix- 
ture of different Bell-states, (ii) Equal weight alternative covariance mag- 
nitude for the Bell-state mixture. The parametrisation in both cases is 
p = x\bi){bi\ + {1 — x)\b2){b2\, where \bi) = \uu ± dd)/^/2, I62) = \ud ± 
du) j The Bell states in these parametrisations may be used in any com- 
bination allowed by the above expression (i.e. all four sign combinations are 
permitted) . The result is not a maximisation over all local unitary transfor- 
mations. 
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p oc 1. This behaviour of the alternative covariance for 'partially entangled' 
configurations is what one would naively expect for a mixed configuration 
entanglement measure E{p). However, there is a particular class of mixed 
states which are disentangled according to the commonly accepted definition 
of a mixed state, whilst exhibiting non-zero alternative covariance properties. 



4.2 Conditions on entanglement measures 

This section assesses the alternative covariance as a measure of entanglement, 
according to the principles outlined by Vedral and Plenio[]I|. 

Firstly, consider the cases where a state exhibits zero quantum correlation. 
A mixed configuration is usually defined to be separable if it can be written 
as a sum of separable projectors; thus 



Then for any such p, 

-tr{(j:P.P^'A^pr^) (Ep^'^pf'B)) 

= Y.Y.p^pMf^'^pf®f^'Bpf 

which is guaranteed to vanish if the p^"'' are projectors. 

Now consider cases where the density matrix is a mixture like 

\\uu){uu\ + \\{u + d){u + d)){{u + d){u + d)\. 
L o 

Here, it can be shown that these mixtures produce non-zero Cp(A, 5), be- 
cause their expansion into projector traces may be used to extract off-diagonal 
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elements of observables A, B in suitable local bases. Take for instance the 
two local operators on separate spin | bases, A = B — oi-, which is a uni- 
tary transformations of the equal- weight matrix, (73 used previously. The 
terms in the expansion of C above are non- vanishing only if non-parallel, 
non-orthogonal projectors are used. Thus we evaluate the two non- vanishing 
terms, with local projectors \u -|- -|- (i|/2, giving 

(it d\u) {u\A\u^ d){{u\B\u + d){u + d\u) -{u + d\B\u){u\u + d)) 

= H)H-0 = -2, 
{u\u + d){u + d\A\u){{u + d\B\u){u\u + d) - {u\B\u + d) {u + d\u)) 

= + = -2. 

These terms do not cancel, and the other two terms in the expansion are 
zero, so the alternative covariance Cp is non-zero even though we are dealing 
with a disentangled state, according to the usual terminology. 

From this we deduce that the measure E based on alternative covariance 
is not a measure of entanglement (according to the conditions provided by 
Vedral and Plenio), when mixed states are encountered, since it violates the 
first condition for such measures. 

4.3 Local purification procedures 

The third condition on entanglement measures proposed by Vedral and Plenio 
is that the entanglement of a state should not increase under the three types 
of purification processes (LGM, CC, and PS). 

Let us consider all possible LGM, CC and PS measurement operations 
which act on a disentangled state 

Po = \uu){uu\ 

and yield states which exhibit non-zero quantum covariance Cp{A^^\ B^'^^), 
such as 

1 1 

pi = -\uu){uu\ + -^\{u + d){u + d)){{u + d){u + d)\ 
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The change from po to pi may be performed by the classically correlated set 
of local measurements Vi . . . Vg, where 



1 1 

Vi — —^\uu){uu\ — ^=\u){u\<^\u){u\ 
V 2 V 2 

V2 = ^\{u + d){u + d)){uu\^^\{u + d)){u\^\{u + d)){u\ 

together with V3, . . . , Vg which have {ud\, {du\, {dd\ in place of {uu\. Thus 
^iPo^i = Pi) J2 ^i^i 0*^ 1) so this represents a complete measurement. 
Because this set of operations Vi is local, we would expect any measure 
of entanglement to not increase under these operations. However, the state 
which results from this procedure docs in fact exhibit non-zero Cp, as was 
demonstrated earlier. Therefore we have found a complete local general 
measurement for which the entanglement increases based upon alternative 
covariance. Thus the unitarily maximised Cp does not satisfy the conditions 
for a measure of mixed state entanglement proposed by Vedral and Plenio. 



5 Conclusion 

We have shown that for pure two 2-statc systems, the maximised covariance 
agrees with other measures of entanglement in specifying which states are 
disentangled and which are maximally entangled. For subspaces of larger 
dimension, the situation is less clear-cut since there are ambiguities in the 
process of selecting eigenvalues for the local operators, but it is possible 
to obtain information about the degree of higher-order correlations of two 
subsystems in this way. 

The variance and covariance could be used to indicate the best measure- 
ments to make to detect entanglement of bipartite systems, by locating the 
unitary transformation which produces maximum covariance. However, the 
problem of mixed state separability measures is not resolved by the corre- 
lation functions C and the covariance, as we have used them in this paper. 
Nevertheless the link between alternative variance and the conjugate vari- 
able derivatives for position, momentum, energy and time is intriguing, and 
might be extended to other quantised models involving conjugate variables. 
The inequality var(X) > C{X,X), with equality for pure states, should be 
applicable to situations where observations are made upon impure states. 
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It is interesting to observe that for pure density matrices p, both the 
average value and the variance of an operator F are symmetrical under the 
interchange of p and F, when expressed in the form 

(F) = tr(pF), var(F) = Cp{F, F) = ii{p^F^ - pFpF). 

They motivate the suggestion that a quantum state \ip) is best represented 
not as a ket, but as a projection operator = \ip){ilj\. Any measurement 
process on the state is relational in the sense of Rovelli the notion of 
state vector reduction is replaced by the notion that every new measure- 
ment process requires a new Hilbert space to be defined, with new operators 
corresponding to the new state and any successive measurement made upon 
the state. The basic idea is that a measurement treats the state informa- 
tion in a comparative sense, with meaning only in relation to the operator 
corresponding to the observable quantity being measured. 

The process of state reduction to the eigenstate of some observable under 
a measurement is an apparent one resulting not from the action of the opera- 
tor corresponding to the measurement, but from some ot/ier physical process 
which occurs in the measurement, such as a filtering process of some kind 
(like the case of a Stern-Gerlach experiment), causing /ater measurements on 
the physical system to be represented in a completely new basis by new state 
and measurement operators. The adoption of the projector as a unit of re- 
lational information does not change any predictions of quantum mechanics 
in terms of average valuesQ 

A The Majorana-Penrose Representation of 
Symmetrised States 

It is well known that a spin state \j, m) can be obtained by forming a fully 
symmetric direct product oi n = 2j spin 1/2 states. Denoting the spin 1/2 

^The paradox of Schroedinger's Cat is resolved by representing the cat's state and the 
measurement on an equal footing as projection operators in the (relational) basis for that 
measurement. Traces are taken to predict average real-number values, but at no stage does 
one say 'the cat is in a state', since the state projection operator only has relational 
significance in this interpretation. 
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states by \u) = |l/2, 1/2), \d) = |l/2, —1/2), the spin j states are given by: 

2j 

|j,j) = l^Iu?^) 

— = [\du . . .u + ud . . .u + ■ ■ ■ + uu . . .d)]/^/2j 
\j,m) = $^|M^'+™(rfy-™)/v/^^ (14) 

sym 

li, -j) = \dd...d). 

In the Majorana-Penrose representation |^ these states are mapped onto the 
unit sphere, with stereographic projection z taken from the South Pole, onto 
the complex plane by making the association: 

|j,m)^ v/^^^^^+™. (15) 

This association may be rewritten in terms of the spin-1 basis in an 
elegant way which emphasises that the 2* powers are proportional to a sym- 
metrised sum where i spin-^ terms of spin |u), and {n — i) spin \d) terms are 
selected, namely 

'^^Cj+mZ^^'^ ^ \U . . .Ud . . .d) + . . . + \d . . .du . . .U) (16) 

where the sum is taken over all ways of selecting i spin-i elements from the 
2j + 1 elements in the spin-j space. 

Thus from a general spin state = ^rn'^jmlJ™) can define a poly- 
nomial 

p{z) = a2jz'^^ + a2j^iz^^^^ H h CqI, 

where Oj+m = \/^^Cj^rn'4'jmZ-'~^"^- The roots of this polynomial^ give 2j 
points on the stereographic plane and, correspondingly, a general spin j state, 
being some superposition over m values, will be described by 2j distinct 
points on the Poincare sphere, obtained by stereographic projection of the 

^The action of exponential functions of angular momentum operators on Majorana- 
Penrose polynomials are amusing. We quote them without proof: (i) exp(^ Jz).p(z) = 
exp(-0')p(^e«); (ii) exp(CJ-).p(z) =p(z + C); (iii) exp(^J+).p(z) = z^^p{z/{Cz + 1)). 
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North pole 




South pole 



Figure 5: Majorana representation of the state with polynomial {z — i){z + i) 

roots placed on the horizontal x-y plane. It is worth remarking that another 
spin state {ip'), represented by another polynomial, 



yields a scalar product, 



H h agl. 



In the case where the degree m of the polynomial is less than 2j, 2j — m 
additional points at the South pole (the projective point) are added to the 
projection, corresponding to 'roots at infinity'. Figure |^ depicts an example 
of stereographic projection for the spin-1 state represented by the quadratic 
polynomial 



p{z) = {z- i){z + i) = - 1 ^ + |1 



2. 



(17) 



By contrast, for spin 1, the maximum and minimum spin states |1,1) and 
1) correpond to two repeated points at the North and South Pole, 
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respectively; on the other hand the pure intermediate spin state |1,0) = 
\ud + du) / \/2 corresponds to one point at the North pole and the other at 
the South pole. The Bell states, |1,±) = \uu ± dd)/\/2 are also antipo- 
dal points, but are located around the equator, specifically at coordinates 
(1,0, 0)& (-1,0,0) and (0,1,0) & (0,-1,0). Under all accepted measures of en- 
tanglement, the spin zero state |0, 0) = \ud — du)/^/2, the Bell states |1, ±) 
and |1, 0) (in fact, all local unitary transformations of these states) are max- 
imally entangled. This suggests that taking the two points as "far apart 
as one another" on the Poincare sphere is one possible way of producing 
maximal entanglement. 

Another point of interest is that for these spin 1 states, the rotationally 
invariant dispersion measure, 

(AJ)2 = -(/).(/), 

attains the maximum value of 2, because (J) vanishes; so it also suggests 
that the quantity (AJ)^/(J)^ might serve0 as a way of characterising the 
entanglement of the individual 1/2 spins that make up the j state. In fact 
the dispersion in J is equivalent to dispersion in any of the components J'-*-' 
and to the covariance of any two components. This is because the total 
angular momentum is J = J^^^ + J^^-* + . . . + J*^"^ and all states on the 
Poincare sphere are symmetrised; therefore the values taken by all the local 
J*^*) are the same. Thus when acting on these symmetrised states, J = nJ^'^'^ , 
for all i = l,2,...,n, and the variance of any local J*^*^ is just as good 
an entanglement measure. We can also take any two subspaces i,j and 
evaluate the covariance of J'-*^ and J'--'^; the space being symmetrised, any i, j 
(including i = j) may be taken! Thus the symmetrised states of maximised 
(AJ)^ have maximal covariance of the local J^^\ as well. Finally, since we 
are dealing with spin 1/2 states, J^*^ have only have two distinct eigenvalues, 
and hence these local operators are directly proportional to the equal weight 
local operators of dimension two defined in section 3. 

We may use these considerations to find the symmetrised states of max- 
imum dispersion for arbitrarily high J - these states will have the maximum 
covariance of any two local J^*\ From the result in section 4, the local covari- 
ance attains the maximum value of 1 whenever the reduced density matrices 
for the subspace on which those local operators jointly act are of maximum 
entropy. For a spin | space, the reduced density matrices must be of the 
form Pi = l[\ui){ui\ + \di){di\]. 
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For general j-values, observe that disentangled states on the Poincare 
sphere are fully factorisable since the states are symmetrised. Clearly a fully 
factorisable state will have zero local operator variance for any operator with 
the local state as an eigenstate. Their Poincare sphere representation simply 
consists of n = 2j repeated roots, since the general symmetrised, factorised 
state is a product of n kets 

{a\u) + b\d))) . . . {a\u) + b\d))) ...^{az + &)" 

because of the 2;-symmetrisation properties for spin n systems. 

Next it is useful to ask which states maximise (A J)^ if this is to serve as a 
possible indication of entanglement. When j = 1, as wc have seen in section 
3, the states are represented on the Poincare sphere by two diametrically 
opposed points, one example being the state \uu + dd)/\/2. When j = 
3/2 the states having maximum (AJ)^ are those states which maximise the 
covariance of any two local operators. As was seen such states must be 
expressible as symmetrised unitary transformations of {\uu) + \dd,))/y^ in 
the local basis for the two operators in question. Thus the overall states must 
be expressible, by symmetry, as global rotations of \uuu+ddd) / \/2. These are 
the 'triangular' states, namely those represented by 3 points on the Poincare 
sphere arranged in an equilateral triangle around any great circle — a global 
unitary transformation (which preserves symmetrisation) simply rotates this 
configuration around the sphere. Choosing the circle to lie equatorially, a 
polynomial producing such roots is p{z) = + 1, which corresponds to the 
state 1-0) = [|3/2, 3/2) + |3/2, -3/2)]/V2 = \ uuu + ddd)/V2. 

When j — 2 the symmetrised states of maximum dispersion are the states 
with a tetrahedral representation on the sphere. Choosing one of the apices 
of the tetrahedron at the North Pole, we arrive at the polynomial 

p(z) = z{z + V2){z- V2e'''^^){z - v^e"^^/^) ^ z^ + 2V2z, 

which corresponds to the maximally entangled state 

1-0) = [|2,2) + V2|2, -l)]/3 = \uuuu + uddd + dudd + ddud + dddu)/V5. 

Note that here, as in the previous case, the maximal lead to vanishing 

— * 

mean values {ip\J\^). 

For j = 5/2 there are two classes of states with maximum dispersion, with 
slightly different geometries. The first of these is pyramidal with one apex at 
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Figure 6: Majorana representation of the state with polynomial p{z) — 



the North pole and the other four apices at equal latitude (or any global rota- 
tion of this state) , while the second has one point at the North pole another at 
the South Pole and the remaining three points distributed equilaterally on the 
equator. Hence the first configuration corresponds to the polynomial p(z) — 
+ z^/V^ or the state [ip) = [|5/2, 5/2) + ^/5/3\5/2, -3/2)]/ ^/sjs, whereas 
the second configuration leads to p(^) = z^ + z, or the state [| 5/2, 5/2) + 
|5/2, — 3/2)]/\/2. Both choices have maximum (AJ)^ and so the total spin 
variance is unable to discriminate between them. However, local equal-weight 
measures are able to discriminate between these states, so it seems that in 
more complicated cases entanglement is most naturally described by keeping 
the covariance of local operators in mind. 

In order to use the (rotationally invariant) dispersion as a measure of 
pure state entanglement, one needs to consider that the symmetrised nature 
of the state is a reflection of the choice of basis. Thus one might define the 
variance-entanglement for a general state as the variance of the symmetrised 
form of the state, under an appropriate local unitary transformation. How- 
ever, it is not always possible to symmetrise an arbitrary state with local 
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unitary transformations in spaces of spin 1 or higher. Thus it seems that the 
dispersion is not a perfect entanglement measure, although it does give an 
indication of the degree of entanglement of these particular states, because 
of its connection with the covariance through symmetrisation. 

B An Integrity Basis for Density Matrix In- 
variants 

Consider p for a composite A^'^^^A^'^^^-dimensional system as an object trans- 
forming under [/(iV^^)) x [/(iV^^)) like {1}{1} x {1}{1}, where denotes 
the reducible N"^ representation of U{N)] in tensor notation we can write p in 
the form p^-^, where early Latin letters refer to the first unitary group and the 
later letters stand for the second group. Here we want to count the number 
of U{N'^^^) X f/(A^(^)) singlets Sn in the symmetrised product p". Thus in the 
notation where representations are labelled by partitions [[T^], 

p" - ({i}{i} X {i}{i}) ® w ^ J2 X i'^}^^}- 

KoXofiou(^{n} 

But {/t}{A} can only contain a singlet if k = A. Moreover k, = X and p = u 
have to be respectively A^*^^^ and A^^^-* part partitions, say "k I~jv(i) n" , etc., 
otherwise {k} vanishes in U{N^^'>). Therefore 

p"|{o}x{o} C J2 {'^H'^} X {A}{A}. 

But it is known that a o a 3 {n} for any a \- n and moreover the order is 
immaterial, because the Clebsch series is symmetric. Thus we have only to 
count the appropriate partitions, 

Sn = |p"l{o}x{o} ^ • ^ ^nw n, A 1-^(2) n}\ . 

This is not easy to work out in the general case, but is relatively simple for 
the case A^*-"^-* = A^^^-* = 2. When n is even or odd, the partitions are: 

n = 2k: {2k,0}, {2k- 1,1},... {k,k} 

n = 2k + l: {2k + 1, 0}, {2k, 1}, . . . , {A; + 1, k}. 
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Now the generating function for invariants of order n in p is written 
-^(5) '■— Yl'^=o^nQ"- Including the even and odd cases, 



The denominator of F{q) is crucial for its interpretation: we can recognize 
that for the 2x2 case invariants are freely generated by two quadratic factors 
(namely {1 — q^Y) and one linear factor (viz. (1 — g)), but there is an extra 
quadratic factor in the numerator which may only be used once. We may 
associate these factors with 

Linear : trp = tr(i)tr(2)p = 1 anyway, 

Quadratic : xi = tr(i) (tr(2)p)^ = p^-p^j , 
& X2 = tr(2) (tr(i)p)^ = 

Extra quadratic : ^aa'^''''' ^w^"' PtjPb'j'- 

The last of these invariants is obviously related to Xi, X2 and tr(p^); however 
we do not get a new invariant from its square because it is then expressible 
entirely as products of Xi)X2, tr and (tr p)^ = 1. Thus effectively tr(p^) 
is only allowed once. 

Wc conclude that any local unitary invariant entanglement measure in 
the 2x2 case must take the form: 

E(p) := F(xi,X2) +tr(p2).G'(xi,X2), 

where F, G are functions which depend on the way E{p) is defined. We think 
that this result must be useful for classifying entanglement. 
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